Analytic treatment of the black-hole bomb 
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A bosonic field impinging on a rotating black hole can be amplified as it scatters off the hole, 
a phenomena known as superradiant scattering. If in addition the field has a non-zero rest mass 
then the mass term effectively works as a mirror, reflecting the scattered wave back towards the 
black hole. In this physical system, known as a black-hole bomb, the wave may bounce back and 
forth between the black hole and some turning point amplifying itself each time. Consequently, the 
massive field grows exponentially over time and is unstable. Former analytical estimations of the 
timescale associated with the instability were restricted to the regimes Mfi 3> 1 and Mfi <C 1, where 
M and fi are the masses of the black hole and the field, respectively. In these two limits the growth 
rate of the field was found to be extremely weak. However, subsequent numerical investigations 
have indicated that the instability is actually greatest in the regime Mfj, = O(l), where the previous 
analytical approximations break down. Thus, a new analytical study of the instability timescale 
' for the case Mfi = O(l) is physically well motivated. In this Letter we study analytically for 

the first time the phenomena of superradiant instability (the black-hole bomb mechanism) in this 
physically interesting regime - the regime of greatest instability. We find an instability growth rate 
of r _1 = m = 1.7 x lO^M" 1 for the fastest growing mode. This instability is four orders of 



on 
o 
o 

(N 
-i— » 

o 
O 



cr 



ON 

o 



magnitude stronger than has been previously estimated. 



Black holes are the fundamental "atoms" of general relativity. They also play a central role in high energy physics, 
astrophysics, and even in condensed matter physics. The fundamental role of black holes makes it highly important 
■ to study the nature of their stability: If a black hole is perturbed in some small way, will the perturbation die away 
over time? Or will it grow exponentially until it can no longer be considered a perturbation and hence demonstrate 
C*~) • the instability of the black hole? 

The issue of black hole stability was first addressed by Regge and Wheeler [l[ who demonstrated the stability of the 
spherically symmetric Schwarzschild black hole. If a Schwarzschild black hole is perturbed, then the perturbation will 
oscillate and damp out over time ■ This implies that perturbation fields would either be radiated away to infinity 
or swallowed by the black hole. 

The stability question of rotating Kerr black holes is a bit more involved. Press and Teukolsky [H, H| have shown that 
rotating black holes are stable under free gravitational perturbations (see also [H and references therein). However, 
. £^ | the superradiance effect may change this conclusion. Superradiant scattering is a well-known phenomena in quantum 
systems [H, Q as well as in classical ones 0, Considering a wave of the form e lm <P e - t ^ t incident upon a rotating 
object whose angular velocity is O, one finds that if the frequency lo of the incident wave satisfies the relation 

lo < mQ , (1) 

then the scattered wave is amplified. 

A bosonic field impinging upon a rotating Kerr black hole can be amplified if the superradiance condition ([1]) is 
satisfied, where in this case f2 = r a ° g2 is the angular velocity of the black-hole horizon. Here r+ and a are the horizon 

radius and the angular momentum per unit mass of the black hole, respectively. The energy radiated away to infinity 
may actually exceed the energy present in the initial perturbation. Feeding back the amplified scattered wave, one 
can gradually extract the rotational energy of the black hole. Press and Teukolsky suggested to use this mechanism 
to build a black-hole bomb (Toj : If one surrounds the black hole by a reflecting mirror, the wave will bounce back and 
forth between the black hole and the mirror amplifying itself each time. Thus, the total energy extracted from the 
black hole will gradually grow. 

Remarkably, nature sometimes provides its own mirror 0]: If one considers a massive scalar field with mass 
M. scattered off a rotating black hole, then for u> < /i = A4G/hc the mass term effectively works as a mirror 
111 Il2l ll3L 14l Il5l| . The physical idea is to consider a wave packet of the massive field in a bound orbit around the 
black hole [12l . 1 1 31 ] . The gravitational force binds the field and keeps it from escaping to infinity. At the event horizon 
some of the field goes down the black hole, and if the frequency of the wave is in the superradiance regime |1| then the 
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field is amplified. In this way the field is amplified at the horizon while being bound away from infinity. Consequently, 
the massive field grows exponentially over time and is unstable [l3j |. 

The nature of the superradiant instability (the black hole bomb) depends on two parameters: 

• The rotation rate a of the black hole. 

• The dimensionless product of the black hole mass M and the field mass ll. The product Mil is actually the 
ratio of the black hole size to the Compton wavelength associated with the rest mass of the field, (we shall 
henceforth use natural units in which G = c = U = 1. In these units fx has the dimensions of 1/length.) 

Former analytical estimates of the instability timescale associated with the dynamics a of massive scalar field in 
the rotating Kerr spacetime were restricted to the regimes Mfx 3> 1 [k| and Mfi <C 1 [13, 14 1. In these two limits 



the growth rate of the field (the imaginary part ioj of the mode's frequency) was found to be very weak, scaling like 



M~ 



-1.84JW/J 



<t 1 for Mil » 1 [12j and like M _1 (M/i) 9 < 1 for the Mfx < 1 case [1J, |l4j. We note, however, that 



the former analytical approximations [l2, 1^, 14| fail in the regime M/i — O(l). Thus, direct numerical integration 
of the perturbation equations seemed necessary to find the actual growth rate of the perturbations in this regime 
Til . 3, EH ■ These numerical investigations have indicated that the superradiant instability is in fact greatest in the 



regime Mfx = O(l) [H, [LH. [l5|. A new analytical study of the superradiant instability in the regime Mfx — 0(1) is 
therefore physically desirable. 

It is worth noting that, previous numerical investigations 12, 14, IH of the black-hole bomb have indicated that 
the instability is most effective under the following conditions: 



• The black hole is maximally rotating with a ~ M. 

• The dimensionless product M/i satisfies the relation M/i 



• The frequency of the unstable mode satisfies the relations to ~ mil and lo ~ fx. (Of course, for the mode to be 
in the superradiant regime one should have lo < mil. In addition, for the mode to be in a bound state it should 
satisfy lo < fx.) 

As we shall show below, the black-hole bomb and the associated instability timescale can be studied analytically 
in the above mentioned regime of physical interest (the regime of the greatest instability) . The physical system we 
consider consists of a massive scalar field coupled to a rotating Kerr black hole. The dynamics of a scalar field W of 
mass \x in the Kerr spacetime 16[ is governed by the Klein-Gordon equation 



One may decompose the field as 



(V Q V a - /i 2 )* = . 



9 lm (t, r, 9, <f>) = e m *Si m (9; auj)R lm (r; aoj) 



(2) 



(3) 



where (t, r, 6, 4>) are the Boyer-Lindquist coordinates [1(J, ui is the (conserved) frequency of the mode, I is the spheroidal 
harmonic index, and m is the azimuthal harmonic index with — I < m < I. (We shall henceforth omit the indices I 
and m for brevity.) With the decomposition (IJ? and S obey radial and angular equations both of confluent Heun 
type coupled by a separation constant A(auj) l(17l Il8j|. The sign of u>i determines whether the solution is decaying 
(u>i < 0) or growing (o>j > 0) in time. 

The angular functions S(0; acu) are the spheroidal harmonics which are solutions of the angular equation 0, [l8| 
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a 2 (u> 2 - fx 2 ) cos 2 
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5 = 0. 



(4) 



The angular functions are required to be regular at the poles 9 = and 9 = n. These boundary conditions pick out 
a discrete set of eigenvalues {A{\ labeled by an integer I. For to ~ /x one can treat a 2 (ui 2 — /x 2 ) cos 2 9 in Eq. ([4]) as a 
perturbation term on the generalized Legendre equation. We can then expand the separation constants in powers of 
a 2 (Li 2 -L0 2 ) to find 



A = l(l + l) + J2c k a 2k (LX 2 -u; 2 ) k . 



k=l 



The expansion coefficients {cfc} are given in [For example, for I 

below), we find c x = -4/5, c 2 = -4/875, c 3 = 8/65625, ...]. 



(5) 



1, the case of physical interest (see 
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The radial Teukolsky equation is given by 

d ( . dR\ 



(a 

dr V dr 



) + 



K 2 - A(a z uj 2 - 2mauj + tfr* + A) 



R = 



(6) 



where A 



2Mr + a 2 and K = (r 2 + a 2 )uj — am. The zeroes of A, r± — M ± (M 2 — a 2 ) 1 / 2 , are the black hole 



(event and inner) horizons. 

We are interested in solutions of the radial equation ([6]) with the physical boundary conditions of purely ingoing 
waves at the black-hole horizon (as measured by a comoving observer) and a decaying (bounded) solution at spatial 
infinity [HQl]. That is, 



ijj^" 2 "" 28 asr^oo (y^oo) ; 

where the "tortoise" radial coordinate y is defined by dy = [(r 2 + a 2 )/A]dr. These boundary conditions single out a 
discrete set of resonances {uj n } which correspond to bound states of the massive field 12. Il5|. 
It is convenient to define new dimensionless variables 

r — r + r + — r_ lo — mf2 

x = ; r = ; H7 = — — ; k = 2ivr+ , (8) 

r+ r + 2-kTbh 



in terms of which the radial equation becomes 



x(x + t)^-4 + (2x + t)^ + VR = , (9) 
oar dx 



where V = K 2 lr\.x{x + r) — \a 2 uj 2 — 2mauj + \j?r 2 Ax + l) 2 + A] and K = r+cox 2 + r + kx + r+mr/2. 



As discussed above, previous numerical investigations [12j, [14|, |15| have indicated that the black-hole instability is 
most pronounced in the regime t«1 with M(mfl - u)<l. As we shall now show, the radial equation is amenable 
to an analytic treatment in this regime of physical interest. 

We first consider the radial equation ([9]) in the far region x 3> max{r, M(mf2 — a;)}. Then Eq. is well 
approximated by 

9 d 2 R dR , x 

x 2 ^j + 2x— + V fal R = 0, (10) 

where Vf ar = (cj 2 - ^ 2 )r\x 2 + 2{u>k - fi 2 r + )r + x - (a 2 uj 2 - 2mauj + [i 2 r\ +A-k 2 ). A solution of Eq. (JTHJ) that satisfies 
the boundary condition |(7J) can be expressed in terms of the confluent hypergeometric functions M(a,b,z) 0, HI] 

R = d(2VM 2 - uo 2 r + )^ +l3 x-^ +t) e-^^- u2r + x M{^ + - k, 1 + 2/3, 2\J fi 2 — u 2 r + x) + C 2 (0 -> -0) , (11) 



where C\ and C 2 are constants and 



2 = a 2 u? - 2mau + u 2 r\ + A - k 2 + ; K = —. '~ = . (12) 

+ 4 ^fjfl _ ^2 

The notation (0 — > —0) means "replace by —0 in the preceding term." 

We next consider the near horizon region x <C 1. The radial equation is given by Eq. ([9|) with V — > V^ear = 
— (a 2 to 2 — 2mauj + y^r 2 ^ + A) + (kx + vjt /2) 2 / x(x + t) . The physical solution obeying the ingoing boundary conditions 
at the horizon is given by 3, 23| 

i /X — fc) I I 

R = x~i TU (- + l) 2 F 1 (-+0-ik, -- 0-ik;l-iw\-x/T) , (13) 

\t I 2 2 

where 2^1(0, b; c; z) is the hypergeometric function. 

The solutions (fTTj) and f| 13|> can be matched in the overlap region max{r, M(mf2 -w)}<x<l. The x <C 1 limit 
of Eq. (HII) yields~0,[23 

i? -» C 1 (2 v / / i 2 -^ 2 r+)H/3 2: -|+/3 + c 2 (0^ -0) . (14) 



1 _ uk — /i 2 r + 



4 



The x > t limit of Eq. (T3J) yields [H, Hj| 



r(2/3)r(i - iw) 



x -2+P + (/3 _> _/3) . 



r(i + /3 - «fe)r(i + /3 - + ifc) 

By matching the two solutions in the overlap region one finds 



Ci =T 



_ i 0-iw/2_ 



r(2/?)r(i - iw) 



r(i + /3 - ik)T{\ + P-izu + ik) 



(15) 



(16) 



C* 2 



,i+/3-iro/2_ 



r(-2/3)r(i 



ITU 



Approximating Eq. (jTTJ) for 



r(i -/3-ifc)r( 
oo one gets 0, S3] 



iw - 



■iife) 



(2^ 



w 2 r+ 



,-4+0 



(17) 



i? 



d(2vV^r + ) 



-k r(i + 2/3) , 
r(± + /?-«)" 



+ C 2 (/3^-/3) 



\/M 2 -^ 2 i , + a: 



■ V VM +; T(i + /3 + K ) 



1+K (-l)-5- /3+K + C , 2 ( / 9^ -/3) 



a --\/ fJ, 2 -u 2 r + x 



(18) 



A bound state is characterized by a decaying field at spatial infinity. The coefficient of the growing exponent 

should therefore vanish. Taking cognizance of Eqs. (fl6|) - ([T8l) for w ^> 1, one finds the 



evV-^'+x in Eq. JlS 
resonance condition for the bound states of the field 



1 



(8iy 



r(-2/3) 



r(2/3) 



r(| + /3-ifr) 



/3- 



ifc)r(i 



/3-«) 



M 2 y/n 2 -u> 2 (mn-uj) 



1 2/3 



(19) 



The growing resonances of the field can be estimated analytically in the regime of physical interest to ~ y ~ mfi: 
We first note that in this regime the right-hand side (RHS) of Eq. (fT9]) is small [due to the factors M^/y 2 — u> 2 
and M(mD, — u))]. One may therefore write a zeroth-order approximation for the resonance condition: LHS = 
l/r(i + P — k) ~ 0. Using the well-known pole structure of the Gamma functions one finds the approximated 
resonance condition 



(20) 



where n > is a non- negative integer. Taking cognizance of Eqs. ([5]) and (fT2|) . one realizes that Eq. (|20|) is a simple 
polynomial equation for the variable M 2 (y 2 — lu 2 ), whose solutions we denote by u>jj.K [Note that for the zeroth-order 



approximation one has ujf^ = 0.] One can then use to)^' in the equation SyLHS^)^, u>\ 
obtain a simple polynomial equation for the first-order solution 

In figure [1] we depict results for the most unstable (fastest growing) mode with I = m = 1. We present results of 
the direct solutions of both the exact resonance condition (fT9l) and the polynomial approximation ([20]) . One finds 



(0) 



,(°) 



to 



a good qualitative agreement between the two. The maximum growth rate we find is r 



1.7 x 10— 



where r is the e-folding time. We would like to emphasize that, this growth rate is four orders of magnitude stronger 
than has been previously found. 

What are the observable consequences of this instability? Let us consider for example the neutral spinless pion 
7T° whose mass is y ~ 134.97 MeV/c 2 . The superradiant instability is most pronounced for My ~ 0.469, which 
corresponds to a primordial black hole of mass M ~ 9.3 x 10 11 Kg [15]. For the superradiant instability to be 
effective, the lifetime of the neutral pion, t\/2 ~ 8.2 x 10~ 17 sec, should be significantly longer than the timescale 
associated with the instability, r = (1.7 x X0~ 3 )~ 1 GM/c 3 ~ 1.3 x 10~ 21 sec. This condition is indeed satisfied by more 
than four orders of magnitude. One therefore concludes that the superradiant instability in the neutral pion channel 
may indeed manifest itself for primordial black holes. 

In summary, we have studied analytically the instability of rotating black holes to perturbations of massive scalar 
fields. Former analytical estimations 13, 13, EH of the timescale associated with the instability were restricted to 
the regimes My 1 and My <C 1. In these two limits the growth rate of the field was found to be extremely weak. 
However, subsequent numerical investigations T3, 14, l3] have indicated that the instability is actually greatest in 
the regime My = O(l), where unfortunately the previous analytical approximations are not suitable to describe the 
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FIG. 1: Superradiant instability for maximally rotating black holes (a ~ M). The results are for the I = m = 1 mode, the mode 
with the greatest instability. The growth rate Mwi is shown as a function of the dimensionless product M[i. We display results 
for the direct solutions of both the exact resonance condition (|19[) and the polynomial approximation (I20|l . The maximum 
growth rate is r" 1 = wi — 1.7 x 10~ 3 M _1 , where r is the e-folding time. 



dynamics of this instability. Motivated by these numerical studies, we have provided here for the first time an analytic 
treatment of the superradiant instability (the black-hole bomb mechanism) in the physically most interesting regime 
M/U ~ |, where the instability is most pronounced. We find an instability growth rate of t _1 = uj = 1.7 x 10~ 3 M _1 
for the fastest growing mode- four orders of magnitude stronger than has been previously estimated. 
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